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In this work we present a new jet model for the non-thermal broadband emission of low-mass mi- 
croquasars. We calculate the contribution of relativistic particles, primary electrons and protons as 
well as secondary muons, charged pions and electron-positron pairs, to the electromagnetic spec- 
trum of the sources. The distribution in energy of all particle species is obtained for an extended, 
inhomogeneous region. We include detailed analysis of particle energy losses, injection, decay 
and escape from the acceleration zone. We also calculate absorption effects due to photon-photon 
annihilation. As an application, we consider the case of XTE Jill 8+480, a well-known low-mass 
X-ray binary in the galactic halo, and we present predictions about its high-energy radiation. 
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1. Introduction 

Microquasars are X-ray binaries that produce jets, whose presence is inferred due to their 
characteristic synchrotron radio spectrum. The detection of non-thermal radiation reveals that mi- 
croquasar jets are able to accelerate particles (at least electrons) up to relativistic energies. Four 
microquasars have been detected in gamma-rays; these very high-energy photons are thought to 
be emitted at the jets. The origin of the gamma-ray emission is believed to be leptonic, but it is 
possible that relativistic protons, or even secondary pions, muons and pairs, may contribute as well. 

In this work we present a jet model for the broadband spectrum of low-mass microquasars 
(LMMQs). These are systems where the companion star is old, dim and without strong winds. 
Since the star cannot provide dense targets, it is expected that in LMMQs any high-energy emission 
is produced by the interaction of the relativistic particles with the magnetic, matter and radiation 
fields of the same jets. 

In our model, acceleration of electrons and protons takes place in an extended, inhomoge- 
neous region of the jets. The energy distribution of relativistic particles is calculated taking into 
account cooling, injection and convection. We then obtain the broadband radiation spectrum due 
to several interaction processes, including the contribution of secondary particles. We also assess 
the relevance of photon absorption through pair creation in photon-photon annihilation. 

Finally, we apply the model to the broadband spectrum of the low-mass X-ray binary XTE 
Jl 118+480. This system is located at a distance ss 1.72 kpc in the galactic halo (galactic latitude 
b « 62.3°), and consists of a K-type star and an 8.5M black hole (Gelino et al. 2006). Although no 
jets have been directly imaged yet, the observed radio emission suggests the presence of outflows 
during the low-hard state (Fender et al. 2001). Applying our model, we fit the observed radio-to- 
X-rays spectrum of this source and make predictions for its high-energy emission. 

2. Jet model 

The model described here is presented in detail in Reynoso & Romero (2009) and Romero & 
Vila (2008). Two symmetrical jets are injected at a distance zo = 50R g from a black hole of mass 
Mbu', Rg = GMbh/c 2 is the gravitational radius of the black hole. The outflow advances with a bulk 
Lorentz factor rj et = 2, expanding laterally as a cone of half-opening angle 0j et « 6°. We assume 
that the jet terminates at z = Zend- The angle between the jet axis and the line of sight is 6 = 30°. 

The accretion power of the black hole is a fraction of its Eddington luminosity, L accr = ^accr^Edd- 
Following the disk-jet coupling hypothesis of Falcke & Biermann (1995), we assume that a fraction 
of the accretion energy powers the jets: 2Lj et = gjet^accr, with q ]tA = 0.1. The factor two takes into 
account the existence of a jet and a counterjet. 

The mean value of the magnetic field at the jet injection point, Bo = B(zo), is determined 
demanding equipartition between the magnetic and kinetic energy densities. For larger z we pa- 
rameterize the magnetic field as a power law, B(z) = Bo(zo/z) 1 ' 5 - 

Part of the jet power is converted into kinetic energy of relativistic particles through a diffusive 
shock acceleration mechanism. The power injected in relativistic particles is L re i = </ re i£jeb with 
q m \ =0.1. Both protons and electrons are accelerated, so L m \ = L e +L p where L p and L e are the 
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power injected in relativistic protons and electrons, respectively. We relate these two quantities as 
L p = aL e , with a > 1 a free parameter. 

The acceleration region is located at z acc > zo- The value of z aC c is such that the magnetic 
energy density is smaller that the jet kinetic energy density at that distance from the black hole (see 
Komissarov et al. 2007 for a discussion on this topic). The injection function of relativistic primary 
protons and electrons is assumed to be of the form 



Here E max (z) is the maximum energy that particles can achieve at fixed z; it is determined by the 
balance of the total energy loss rate and the acceleration rate (see below). The function f(z) is 
given by a step-like function, that is unity for z < Zmax and becomes negligible for larger z- 

Neutral and charged pions are created through the interaction of relativistic protons with 
photons and non-relativistic protons. Whereas neutral pions decay into photons, charged pions 
yield neutrinos, muons and electrons/positrons; muons in turn decay into neutrinos and elec- 
trons/positrons. Electron/positron pairs are also directly injected in proton-photon collisions. For a 
complete list of expressions regarding the injection functions of secondary particles, see Reynoso 
& Romero (2009), Romero & Vila (2008) and references therein. 

The isotropic steady-state energy distribution of relativistic primary and secondary particles in 
the jet comoving reference frame, N(E,z) (erg -1 cm -3 ), is calculated solving the transport equation 
(Khangulyan et al. 2008) 



This equation takes into account particle injection, convection, energy losses and removal of parti- 
cles due to decay. The convective term v conv dN/dz, where v conv ~ vj e t is the particle convection ve- 
locity, allows to incorporate the effect of the variation with z of the parameters that characterize the 
acceleration region. This is not accounted for in one-zone models where the acceleration/emission 
region is assumed to be homogeneous, or at least thin enough to neglect all spatial dependence. 
The decay term is non-zero only for pions and muons. 

We include both adiabatic and radiative energy losses. In the case of leptons (primary electrons, 
secondary electron/positron pairs and muons) we consider losses due to synchrotron radiation, rela- 
tivistic Bremsstrahlung and inverse Compton scattering. For protons and pions we take into account 
cooling due to synchrotron radiation and inelastic collisions with photons and thermal protons. The 
synchrotron field of primary electrons is used as target photon field for inverse Compton scattering, 
proton-photon and pion-photon collisions. Convenient expressions for the energy loss rates for all 
these processes can be found in Reynoso & Romero (2009) and Romero & Vila (2008). 

The maximum energy of primary protons and electrons at fixed z is determined equating the 
total energy loss rate with the acceleration rate dE/dt\ acc = T]ecB(z). The acceleration efficiency 
r\ < 1 depends on the detailed physics of the acceleration mechanism; we simply regard it as a free 
parameter of the model. 




[Q] = erg 1 cm 3 s 



(2.1) 




(2.2) 
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For a given set of values of the model parameters, we calculate the radiative output due to all 
the interaction processes mentioned above. The correspondent formulas for the photon emissivities 
and other details about the calculations are presented in Romero & Vila (2008), Reynoso & Romero 
(2009) and references therein. 

Finally, we estimate the correction to the primary emission spectrum due to internal photon 
absorption through pair creation, 77 — > e + e~ . The opacity T rr for this process is calculated as in 
Gould & Schreder (1966). 

3. Results 

Figures [I] to |5] show the results obtained for two generic sets of values of the model parameters, 
cases A and B in Table [l[ The cooling rates for primary protons and electrons in model A are plotted 
in Figures [l] and ||. Synchrotron and adiabatic losses are the main cooling mechanisms of electrons 
for all z, whereas adiabatic losses always dominate for protons. Cooling due to those processes 
that involve interaction with the synchrotron photon field becomes negligible for z > Zmax- This is 
because in this region the number of high-energy electrons decreases drastically (see below) and 
consequently the synchrotron emissivity is quenched. 

The evolution of the particle distributions for model A is shown in Figure ^} At fixed z the 
distributions depend on energy as a power law: N e oc E~ 3 for electrons and N p °< E p 2 in the case 
of protons. For z > Zmax, where Q m 0, there is no acceleration and particles only lose energy. The 




Figure 1: Acceleration and cooling rates for primary relativistic electrons, calculated for the parameters of 
model A at different distances from the black hole. From left to right: z acc , z max and z e nd- 




Log ]n (E/ e V) Log^E/eVl L Dg|n (E/ e V) 



Figure 2: Idem Figure [j] for relativistic primary protons. 
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Table 1: Values of the relevant parameters of the model. 



Parameter 


Model A 


Model B 


XTE 2000 


XTE 2005 




1 Al J 

1\JMq 


1 A A/f 
IUMq 


O.JMq 


O.DM0 


^acc 


1.6 x 10 9 cm 


1.6 x 10 9 cm 


3 x 10 8 cm 


6 x 10 8 cm 


Zmax 


10 11 cm 


10 11 cm 


3 x 10 9 cm 


6 x 10 9 cm 


£end 


10 13 cm 


10 13 cm 


10 12 cm 


10 12 cm 


<7accr 


0.1 


0.1 


0.2 


0.1 




0.2 


io- 2 


10~ 2 


io- 2 


a 


100 


1 


1 


5 


a 


2 


2 


1.5 


1.5 


yM 
'mm 


10 


10 


50 


50 



(*' Minimum Lorentz factor of primary protons and electrons. 




Figure 3: Distributions of relativistic primary electrons (left) and protons (right) as a function of energy and 
distance to the black hole, for the parameters of model A. 



number of high-energy electrons decreases since these cool much faster than the low-energy ones. 
The same is true for protons, but the effect is less pronounced since proton cooling is slower. 

A number of spectral shapes can be obtained varying the value of the parameters. Figure [| 
shows the spectral energy distributions (SEDs) for the parameters of models A and B. In case B, 
where a = 1 and v\ = 10~ 2 , synchrotron and IC radiation from primary electrons completely dom- 
inate the spectrum, reaching a maximum luminosity L y m 10 36 erg s _1 . The high-energy emission 
due to pp interactions is negligible. Case A corresponds to a = 100 and v\ =0.1. The power 
injected in primary electrons is 100 times smaller than in the previous model, and then the syn- 
chrotron luminosity reaches only Ly rj 10 34 erg s _1 . This in turn reduces the IC emission levels 
since there are considerably less target photons. For E y > 1 GeV the dominant contribution is due 
to pp interaction, with a luminosity Ly ~ 10 erg s _1 . Radiation from secondary particles is not 
significant in any case. 

The attenuation factor exp(— T ry ) as a function of photon energy and z is shown in Figure 
|5[ Absorption is only important for photons with Ey > 1 TeV, that annihilate against low-energy 
photons of the electron synchrotron field. The effect is therefore stronger in model B. The effect 
on the SED is, however, not significant in any case. 
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Figure 4: Spectral energy distributions obtained for two different sets of model parameters. In the syn- 
chrotron spectra, e~ , e (77) and e ± (py) indicate primary electrons and pairs injected through photon-photon 
annihilation and direct pair creation in proton-photon collisions, respectively. 




Figure 5: Attenuation factor exp(— lyy) for models A (left) and B (right). 



We applied the model to fit the observed spectrum of the LMMQ XTE Jl 1 18+480 in the low- 
hard state during the outbursts of 2000 (McClintock et al. 2001) and 2005 (Zurita et al. 2006). Here 
we also included the emission of a standard thin accretion disk. The best fit models were obtained 
for the parameters listed in Table [TJ the correspondent SEDs are shown in Figure ^. According to 
these results, the source might be marginally detected by Fermi in the GeV range during a low- 
luminosity low-hard state such as the case of the 2005 outburst. However, during a more luminous 
low-hard state as in 2000, also the hadronic contribution to the emission spectrum at TeV energies 
might be detected with MAGIC II and the future array CTA. 

4. Conclusions 

We have developed a model for the broadband spectrum of jets in LMMQs that includes the 
radiative contribution of electrons, protons and secondary particles generated by several interaction 
processes. The energy distributions of all species were obtained solving the transport equation 
in an inhomogeneous region, taking into account cooling, injection, convection and decay. We 
also assessed the effect of photon absorption in the jet. We applied the model to fit the observed 
spectrum of the low-mass X-ray binary XTE Jill 8+480 and make predictions for its high-energy 
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Figure 6: Fits to the observed spectrum of XTE Jl 1 18+480 in low-hard state during the outbursts of 2000 
(left) and 2005 (right). 

emission. According to our results, luminous low-hard state outbursts of this source might be 
detected at GeV and TeV energies with current and future gamma-ray instruments. 
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